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Abstract. The nature of fast magnetoacoustic and Alfven waves is investigated in a zero /3 plasma in the neighbourhood of a pair 
of two-dimensional null points. This gives an indication of wave propagation in the low j3 solar corona, for a more complicated 
magnetic configuration than that looked at by McLaughlin & Hood (2004). It is found that the fast wave is attracted to the 
null points and that the front of the wave slows down as it approaches the null point pair, with the wave splitting and part of 
the wave accumulating at one null and the rest at the other. Current density will then accumulate at these points and ohmic 
dissipation will then extract the energy in the wave at these points. This suggests locations where wave heating will occur in the 
corona. The Alfven wave behaves in a diff'erent manner in that the wave accumulates along the separatrices. Hence, the current 
density will accumulate at this part of the topology and this is where wave heating will occur. However, the phenomenon of 
wave accumulation at a specific place is a feature of both wave types, and illustrates the importance of studying the topology of 
the corona when considering MHD wave propagation. 

Key words. MHD - Waves - Sun: corona 



1. Introduction 

The solar corona is dominated by the magnetic field. In or- 
der to understand the myriad of phenomena that occur on the 
Sun, it is important to understand the structure (topology) of 
the magnetic field itself. Potential field extrapolations of the 
coronal magnetic field can be made from photospheric mag- 
netograms. Such extrapolations show the existence of an im- 
portant feature of the topology; null points. Null points are 
points in the field where the Alfven speed is zero. Detailed in- 
vestigations of the coronal magnetic field, using such potential 



. field calculations, can be found in |Beveridge et al. (2002)| and 
■ [Brown & Priest (200 1)[ 

. Another key feature of the magnetic topology is the sepa- 
' rator\ a special field line that connects two null points. A 
stable magnetic separator is the intersection of two separa- 
trix surfaces. These surfaces divide the magnetic field into re- 
gions of diflTerent connectivity. A state is topologically stable if 
any small perturbation in the parameters does not change the 
state dBrown & Priest, 1999a| ). Separators have an important 
role in magnetic reconnection ( [Priest & Titov, 199"6| ) and cur- 
rent sheets can form along them (ILongcope & Cowley, 1996). 
McLaughlin & Hood (2004) found that for a single 2D null 
point, the fast magnetoacoustic wave was attracted to the null 
and the wave energy accumulated there. In addition, they found 
that the Alfven wave energy accumulated along the separatrices 
of the topology. The aim of this paper is to see if their ideas 
carry through to magnetic configurations of more than one null 



point. It can be argued that null points appear in pairs, (e.g. 
by a local bifurcation) which makes it relevant to investigate 
multiple null point topologies. A double null point may arise as 
a bifurcation of a single 2D null point ( [Galsgaard et al, 1996 , 
[Brown & Priest, 1998) . 

Waves in the neighbourhood of a single 2D null point have been 
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investigated by various authors. [Bulanov & Syrovatskii (1980) 
provided a detailed discussion of the propagation of fast and 
Alfven waves using cylindrical symmetry. In their paper, har- 
monic fast waves are generated and these propagate towards the 
null point. However, the assumed cylindrical symmetry means 
that the disturbances can only propagate either towards or away 
from the null point. [Craig & Watson (1992)[ mainly consider 
the radial propagation of the m = mode (where m is the az- 
imuthal wavenumber) using a mixture of analytical and numer- 
ical solutions. In their investigation, the outer radial boundary is 
held fixed so that any outgoing waves will be reflected back to- 
wards the null point. This means that all the energy in the wave 
motions is contained within a fixed region. They show that the 
propagation of the m = wave towards the null point generates 
an exponentially large increase in the current density and that 
magnetic resistivity dissipates this current in a time related to 
log 7]. Their initial disturbance is given as a function of radius. 
In this paper, we are interested in generating the disturbances 
at the boundary rather than internally. Craig and McClymont 
(1991, 1993) investigate the normal mode solutions for both 
m = and m ^ modes with resistivity included. Again they 
emphasise that the current builds up as the inverse square of 
the radial distance from the null point. However, attention was 
restricted to a circular reflecting boundary. 
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Multiple null point topologies have also been investigated by 
various authors. [Galsgaard et al, 1996 ) looked at the dynamic 
reconnection properties of a 3D double null. They investigated 
responses of the magnetic field to specific perturbations on the 
boundaries of a 3D box; they tilted the spine axis to form a 
pure m = 1 mode. They found that for a perturbation parallel 
or orthogonal to the separator (between the two nulls), current 
accumulation occured in the separator plane perpendicular to 
the direction in which the spine was moved. For any other ori- 
entation, the current was focussed along the separator line. In 
the model, the boundary motions move the field lines but do 
not return them to their original positions. Thus, the Poynting 
flux induced by the imposed motion (and then fixing the field 
after the motion is complete) accumulates at the resulting cur- 
rent sheet and provides the energy released in the reconnec- 
tion. However, if the boundary motions are simply due to the 
passing of incoming waves through the boundary, then it is not 
clear that the current sheet will form. If this is the case, will 
the energy in the wave, again due to the Poynting flux through 
the boundary, dissipate or simply pass through one of the other 
boundaries? 

Galsgaard & Nordlund (1997) found that when magnetic struc- 
tures containing many null points are perturbed, current ac- 
cumulates along separator lines. The perturbations used were 
random shear motions on two opposite boundaries. Again, 
the field lines were not returned to their original position. 
[Galsgaard et al. (1997a)| looked at shearing a 3D potential 
null point pair, with continuous (opposite) shear on two op- 
posite boundaries (parallel to the separator). This generated 
a wave pulse that travelled towards the interior of the do- 
main from both directions, and resulted in current accumu- 
lation along the separator line with maximum value at the 
null points. From these experiments (which share a lot with 
[Galsgaard et al, 1996} it was concluded that to drive current 
along the full length of the separator line, the perturbation 
wavelength had to be longer than the length of the separa- 
tor line. [Galsgaard et al. (1997b)] looked at perturbations in 3D 
magnetic configurations containing a double null point pair 
connected by a separator. The boundary motions used were 
very similar to those described above (i.e. shear the boundary 
and fix). Their experiments showed that the nulls can either 
accumulate current individually or act together to produce a 
singular current collapse along the separator. However, in all 
these previous works the boundary conditions used have tried 
to mimic the eff'ect of photspheric footpoint motions by moving 
the boundary and holding it fixed. This paper will look at the 
2D null point pair and investigate the propagation and transient 
behaviour of an individual wave pulse entering the magnetic 
structure. 

The propagation of fast magnetoacoustic waves in an in- 
homogeneous coronal plasma has been investigated by 
[Nakariakov & Roberts (1995) , who showed how the waves are 
refracted into regions of low Alfven speed. In the case of null 
points, it is the aim of this paper to see how this refraction pro- 
ceeds when the Alfven speed actually drops to zero. 
The paper has the following outline. In Section [21 the basic 
equations are described. The results for an uncoupled fast mag- 
netoacoustic wave are presented in Section [3 This section dis- 



cusses fast wave propagation with a pulse coming in from the 
top and side boundaries, and numerical and anaytical results 
are presented. Section [4] discusses the propagation of Alfven 
waves and the conclusions are given in Section [5] 



2. Basic Equations 

The usual MHD equations for an ideal, zero jS plasma appro- 
priate to the solar corona are used. Hence, 



p(g.(vV)v 
~dt 



= -(VxB)xB 
= V X (v X B) , 



^+V.(pv) = 0, 
ot 



(1) 
(2) 
(3) 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field), jj. = 
An X 10"^Hm"^ the magnetic permeability, and cr the electrical 
conductivity. The gas pressure and the adiabatic energy equa- 
tion are neglected in the zero p approximation. The magnetic 
diff'usivity is neglected in the ideal approximation. 

2.1. Basic equilibrium 

The basic magnetic structure is taken as a simple 2D two null 
points configuration. There are two such configurations to con- 
sider; one containing a separator and one that does not. The aim 
of studying waves in a 2D configuration is one of simplicity; 
the individual eff'ects are much easier to identify when there is 
no coupling between the fast and Alfven modes. The magnetic 
field with a separator is taken as: 



^0 /„2 



Bo = ^(. 



■ a^, 0, -2xz) , 



(4) 



where is a characteristic field strength and a is the length 
scale for magnetic field variations. This configuration can be 
seen in Figure [T] (left). The other magnetic configuration con- 
sidered takes the form: 



Bo = ^{2xz,0,x^-z^-a^), 



(5) 



and can be seen in Figure [T] (right). Obviously, these magnetic 
configurations are no longer valid far away from the null points, 
as the field strength tends to infinity. 

As in [McLaughlin & Hood (2004)| the linearised MHD equa- 
tions are used to study the nature of wave propagation near the 
two null points. Using subscripts of for equilibrium quanti- 
ties and 1 for perturbed quantities, the linearised equation of 
motion becomes 



the linearised induction equation 

— ^ = V X (vi X Bo) , 
at 



(6) 



(7) 
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Fig. 1. (left) Magnetic field configuration containing two nulls joined by a separator, (right) Magnetic field configuration contain- 
ing two nulls not joined by a separator. 



and the linearised equation of mass continuity 



^+V.(poVi) = 0. 
ot 



(8) 



We will not discuss equation ([8]) further as it can be solved 
once we know Vi. In fact, it has no influence on the mo- 
mentum equation (in the zero p approximation) and so in ef- 
fect the plasma is arbitrarily compressible (Craig & Watson, 
1992). We assume the background gas density is uniform and 
label it as po. A spatial variation in po can cause phase mixing 



( Heyvaerts & Priest, 1983). 

We now consider a change of scale to non-dimensionalise; let 
vi = vv;. Bo = ^oBJ, Bi = 5oBp x = ax\ z = az\ V = ^V* 
and t = tf, where we let * denote a dimensionless quantity and 
y, Bq, a and r are constants with the dimensions of the variable 
they are scaling. We then set = v and \ = a/t (this sets v as 



a sort of constant background Alfven speed). This process non- 
dimensionalises equations Q and (|7]), and under these scalings, 
f = I (for example) refers to ^ = /" = a/v; i.e. the time taken 
to travel a distance a at the background Alfven speed. For the 
rest of this paper, we drop the star indices; the fact that they are 
now non-dimensionalised is understood. 

The ideal linearised MHD equations naturally decouple into 
two equations for the fast MHD wave and the Alfven wave. 
The slow MHD wave is absent in this limit and there is no 
velocity component along the background magnetic field (as 
can be seen by taking the scalar product of equation © with 
Bo. 

The linearised equations for the fast magnetoacoustic wave are: 



dt 



dV_ 



dt 



dV_ 

dx 



(9) 



where the Alfven speed, va {x, z), is equal to + = 



■2(x2-z2) + l 



B, 



(bx, 0, b^) and the variable 



V is related to the perpendicular velocity; V = |^(vi x Bo) • e^j. 
These equations can be combined to form a single wave equa- 
tion: 



dfi 



(10) 



The linearised equations for the Alfven wave, with vi = 
(O, Vy, O) and Bi = (o, by, o) are: 

dvy dby dby dby dvy dvy 
^ = B^c^ B,—^ , —^=B^-^+B,—^, (11) 
dt ox oz ot ox oz 

which can be combined to form a single wave equation: 
d\ Id 

^ = K^^^^)^- '''' 

It is worth noting that the zero /3 approximation will be invalid 
near the null points, where the pressure terms become impor- 
tant. In fact, a finite /3 model would introduce slow waves, 
which would be driven by the excited fast waves. Also, the 
model is linear but the inhomogenity of the medium would lead 
to nonlinear coupling of the modes ( |Nakariakov a/. , 1997| ). 

3. Fast wave 

We consider the eff'ect of sending a fast magnetoacoustic wave 
into these two magnetic configurations from the top and side 
boundaries. Unlike [McLaughlin & Hood (20 04), we also con- 
sider the eff'ect of sending in a wave pulse from a side boundary 
because the magnetic configuration is not symmetric. In fact, 
since the fast magnetoacoustic wave can cross fieldlines, it be- 
haves identically across both magnetic configurations (Figure 
[T] left and right) and so we only have two cases to investigate, 
i.e. ([Tot is true for both magnetic configurations (unlike ([T2l) . 
which depends upon the form of B). 
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3.1. Upper boundary 

We solve the linearised MHD equations ([9]) for the fast wave 
numerically using a two-step Lax-Wendroff scheme. The nu- 
merical scheme is run in a box with -3 < x < 3 and -4 < z < 2, 
with our attention focused on -2 < x < 2 and -2 < z<2. For a 
single wave pulse coming in from the top boundary, the bound- 
ary conditions are taken as: 



V(x, 2) 



sin ojt for <t < - 

CO 

otherwise 



dv 

dx 



dx 



-<>■ T 

oz 



= 



Tests show that the central behaviour is largely unaffected by 
these choices of side and bottom boundary conditions. The 
other boundary conditions on the perturbed magnetic field fol- 
low from the remaining equations and the solenodial condition, 
V-Bi =0. 

We find that the linear, fast magnetoacoustic wave travels to- 
wards the neighbourhood of the two null points and begins to 
wrap around them. Since the Alfven speed, va (x, z), is spa- 
tially varying, diff'erent parts of the wave travel at diff'erent 
speeds, and it travels faster the further it is away from the 
null points (va(x,z) = at the null points). This is a sim- 
ilar eff'ect to that described it the case of a single null point 



in [McLaughlin & Hood (2QQ4)[ However, in this case there is 
a non-zero Alfven speed between the two null points, i.e. for 
-l<x<l,z = 0. The fast wave slows down greatly when 
it approaches this area, but manages to cross this line (the sep- 
arator if we were considering a setup such as that of Figure [T] 
(left)). Hence the wave bends between the two null points and 
passes through the area between them. This part of the wave 
is then again suceptible to the refraction eff'ect, and so contin- 
ues to wrap around the null points, breaking into two waves 
along the line x = (due to symmetry). Each part of the wave 
then continues to wrap around its closest null point, repeatedly, 
eventually accumulating at the null points (x = - 1 and x = -h 1 , 
along z = 0). This can be seen in the shaded contours of Figure 

El 



3.2. Analytical results 

We can approximately solve equation ([TOb for the fast wave to 
gain more insight into the numerical simulations. Substituting 
V = ^''^^'^'^^ • e''"^^ into ^ and assuming that oj » 1 
(WKB approximation), leads to a first order PDF of the form 



— ^ = 0. Applying the method of characteristics. 



dx'' dz 



we generate the equations: 
d(f) 



ds 
dp 

ds 
dq 

ds 
dx 

ds 
dz 
ds 



2x{x^ +Z^ -\)[p^ +q^) 
2z(jc2+z' + l)(/ + ^2^ 
= -P (^'+^'f -2(x'-z') + l 
^ = -gUjc'^z'f-2(x'-z')^l 



where p = q = oj is the frequency of our wave and s is 
some parameter along the characteristic. 
These five ODFs were solved numerically using a fourth-order 
Runge-Kutta method. Contours of constant can be thought 
of as defining the positions of the edges of the wave pulse, 
i.e. with correct choices of s, the WKB solution represents the 
front, middle and back edges of the wave, s is comparable to t 
and so the numerical and analytical work can be directly com- 
pared. The agreement between the analytic model and the lead- 
ing edge of the wavefront is very good, as seen in an overplot 
of a numerical simulation (shaded area) and our WKB solution 
(thick lines) in Figure [2l 

We can also use our WKB solution to plot the particle paths of 
individual elements from the initial wave. In Figure [3l we see 
the spiral evolution of elements that begin at points x = -2, 
-L5, -1, 0, 0.3 0.5 along z = 2. Note that an element that 
begins at x = 0, z = 2 is not aff'ected by the null points and 
passes between the two and off' to -oo. In eff'ect, the two null 
points cancel out each other's eff'ect at this point, although any 
perturbation to an element travelling along this path would send 
the element spiralling towards one of the nulls. The behaviour 
of such a particle can be shown: 



Va(0,Z): 



Hence, solving for z, we find z = tan (A - t), where A = 
arctanzo and zo is some starting position (which is zo = 2 in 
our simulations). 




-2-10 1 2 

X 

Fig. 3. Plots of WKB solution for a wave sent in from the upper 
boundary and its resultant particle paths for starting points of 
X = -2, -1, 0, 0.3 and 0.5 along z = 2. 



3.3. Side boundary 

We also investigated the eff'ect of solving the linearised MHD 
equations for a fast wave coming in from the side bound- 
ary, along X = -2. Again a two-step Lax-Wendroff' numerical 
scheme was run in a box with -2 < x < 4 and -3 < z < 3, with 
our attention again focused on -2 < x < 2 and -2 < z < 2. 
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Fig. 2. Comparison of numerical simulation (shaded area) and analytical solution V for a Fast wave sent in from upper boundary 
for -3 < X < 3, and its resultant propagation at times {a) t=025, (b) t=0.5, (c) t=0J5, (d) t =1.0, (e) t=l.5 and (/) ^=1.75, (g) 
t=2.0, (h) t=2.25, (0 t=2.5, (j) t =3.0, (k) t=3.5 and (/) t=3.15 , labelling from top left to bottom right. The lines represent the 
front, middle and back edges of the WKB wave solution, where the pulse enters from the top of the box. 



For a single wave pulse from the side boundary, the boundary 
conditions were: 



V(-2,z) 



sin cot for <t < - 

CO 

otherwise 



dv 

dx 



x=4 



dV 



0, 



z=-3 



dV 



z=3 



Again, tests show that the central behaviour is largely unaf- 
fected by these choices of the remaining boundary conditions. 

We find that the linear, fast wave travels in from the left hand 
side of the box and begins to feel the eff'ect of the left hand 
side null point (at x = -1, z = 0). The wave thins and begins 
to wrap around this null point. As the ends of the wave wrap 
around behind the left null point, they then become influenced 
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Fig. 4. Plots of WKB solution for a wave sent in from the side boundary and its resultant particle paths. Top left displays starting 
points of z = -2, 0, 0.5 and +1 along x = -2. Bottom left displays a starting point of z = +1.5. The top right subfigure shows a 
starting point of z = 1.86, the lower right shows z = V3. Note the change of scale in the two right hand subfigures. 



by the second, right hand side null point (at x = +1, z = 0). 
These arms of the wave then proceed to wrap around the right 
null point, flattening the wave. Furthermore, the two parts of the 
wave now travelling through the area between the null points 
have non-zero Alfven speed, and so pop through (the wave in 
the lower half plan now travels up, crossing the line -1 < x < 1 
and z = 0, and visa versa for the wave in the upper half plan). 
These parts of the wave break along x = and then proceed 
to wrap around the null point closest to them. Eventually, the 
wave accumulates at both null points. This can be seen in the 
shaded contours of Figure [H 

3.4. Analytical results 

As before, equations [13] were solved numerically using a 
fourth-order Runge-Kutta method, but this time using diff'er- 
ent initial conditions. The thick black lines in Figure [5] shows 
constant cp at difl'erent values of the parameter s. Constant cp can 
be thought of as defining the position of the edge of the wave 
pulse, i.e. with correct choices of s, the WKB solution repre- 
sents the front, middle and back edges of the wave. As before, 
the agreement between the analytic model and the wavefront is 
very good and can be seen in an overplot of a numerical simu- 
lation and our WKB solution in Figure [51 
We can also use our WKB solution to plot the particle paths of 
individual elements from the initial wave. In Figure [H we see 
the spiral evolution of elements that begin at points z = -2, 0, 
0.5, 1.0, 1.5 and 1.86 along x = -2. Note that for a wavefront 



extending between -2 < z < +2, there is a critical distance 
where some elements of the wave are sucked into the nearest 
(left hand) null point, and others spiral round it at such a great 
distance that they feel the influence of the second (right hand) 
null point and begin spiralling towards it (x = +1, z = 0). Note 
that for a starting position of z = V3, the wave element orbits 
in a circle (radius V3). This is the critical starting distance; a 
starting distance, say zo, less than this will sprial into the left 
hand null point, whereas the wave particle will spiral into the 
right hand null point if zo > V3. This can be seen in Figure [H 
Simulations show that this critical distance is - 1, where 
L is the distance between the origin and the side boundary (here 
L = \xo\ = 2). When xo = -2 and zo = V3 (i.e. critical starting 
distance), the equations are satisfied by: 



-1 



z = 



(2+ VSsin^) 
VScos^ 



V3 



cos^ , 



271 . 

, q = TT -\- — — sm 6 , 
(2+V3sin6/) ^f3 



where 6 = 47rs - 7r/3, and confirm the periodic circular orbit. 
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Fig. 5. Comparison of numerical simulation and analytical solution of V for a Fast wave sent in from left hand boundary for 
-3 < X < 3, and its resultant propagation at times {a) t=0.25, (b) t=0.5, (c) t=0J5, (d) t =1.0, (e) t=1.5 and (/) ^=1.75, (g) t=2.0, 
(h) t=2.25, (0 t=2.5, (j) t =3.0, (k) t=3.5 and (/) t=3.15, labelling from top left to bottom right. The lines represent the front, 
middle and back edges of the wave, where the pulse enters from the side of the box. 
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4. Alfven Wave 

The equations describing the behaviour of the Alfven Wave, 
equations (fTTI) . were solved numerically using a two-step Lax- 
Wendroff scheme. [McLaughlin & Hood (2QQ4)| showed that in 
the neighbourhood of a two-dimensional X-point, the linear 
Alfven wave spreads out along the field lines, eventually accu- 
mulating along the separatrices. Hence, for our two null points, 
we hsiwe four cases to consider; an Alfven wave pulse coming 
in from the top boundary and coming in from the side boundary 
for both magnetic configurations shown in Figure [T] 

4.1. Two null points connected by a separator 

We consider the two null point magnetic configuration contain- 
ing a separator, as shown in Figure [T] (left). 

4.1.1. top boundary 

We initially consider a box (-2 < x < 2, -2 < z < 2) with 
a single wave pulse coming in across part of the top boundary 
(I < X < 2, z = 2). We choose such a pulse because, as initial 
experiments showed, the Alfven wave spreads out along the 
field lines as it propagates and we found that this choice of 
boundary condition illustrated this eff'ect much clearer. In fact, 
the wave crosses the separator (z = 0) and so the final, high 
resolution run was performed in a box (-2.5 < x < 2.5, -0.5 < 
z <2). The full boundary conditions were; 



Vy(x, 2) = 2 (sin cot) (sin 7ix/2) for 




dVy 

dx 



= 0, 



dvA 

dVy 

dx 



\<x<2 
otherwise 



= 0, 



-2.5 



dVy 



= 0. 



Tests show that the central behaviour is unaff'ected by these 
choices. The other boundary conditions follow from the re- 
maining equations and the solenodial condition, V • Bi =0. 
We found that the linear Alfven wave travels down from the top 
boundary and begins to spread out, following the field lines. As 
the wave approaches the separator along z = and the sepa- 
ratrix that passes through x = - 1 , the wave thins but keeps its 
original amplitude. The wave eventually accumulates very near 
the separator and separatrix. This can be seen in Figure [6l 

4.1.2. analytical results 

As in all four cases, the Alfven equations we have to solve are: 



d d 
= \B,— ^B,- 
ox oz 



(13) 



In this experiment, = x^ - z - I, = -2xz and we can 
apply the WKB approximation. As in Section [T2l substituting 
= . e~^^^ into ([T3]) and assuming that cl) » 1 , leads to 

a first order PDF of the form (x, z, 0, ^, ^) = 0. Applying 
the method of characteristics, we generate the equations: 

J0 2 

ds 



dp 
ds 



— = (2zq-2xp)^ 



^ = (2zp-\-2xq)^ 
ds 



(14) 



where ^ = [(x^ - z^ - l^p - 2xzqj, p = q = oj is the 
frequency of our wave and s is some parameter along the char- 
acteristic. 




-2 -1 





X 
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Fig. 7. Plots of WKB solution for a wave sent in from the upper 
boundary and its resultant particle paths (thick lines). Starting 
points of x=l, 1.25, 1.5, V7/3 and 1.75 are plotted. 



These five ODFs were solved numerically using a fourth-order 
Runge-Kutta method. Constant can be thought of as defin- 
ing the position of the edge of the wave pulse, so with correct 
choices of s the solution can be directly compared to our nu- 
merical solution (seen in Figure [6] as the dark overplot lines). 
Figures [51 [H [9] and [TOl all show comparisons of the numerical 
simulation and corresponding WKB solutions. In each case, the 
agreement between the analytic model and numerical results is 
very good. 

We can also use our WKB solution to plot the particle paths 
of individual elements from the initial wave. In Figure [71 start- 
ing points of x=l, 1.25, 1.5, ^|1 1?> {separatrix) and 1.75 are 
plotted. 

4.1.3. side boundary 

We now consider a box (-2 < x < 2, -2 < z < 2) with a single 
wave pulse coming in across part of the side boundary (x = -2, 
-1 <z< 1). In fact, the wave never passes the (left) separatrix 
and so the final, high resolution run was performed in a box 
(-2 < X < 0.5, -2.5 < z < 2.5). The full boundary conditions 
were; 



, ^ , f 2(sina;0fsin^l for 
Vy(-2,z) = \ L 2 J 

[ oth 



-1 <z< 1 
< r < ^ 

OJ 

otherwise 
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Fig. 6. Comparison of numerical simulation (shaded wave) and analytical (WKB) solution at times (a) ^=0.25, (b) t=0.15. 
t=l.25, (d) t=lJ5, (e) t=2.25 and (/) t=2J5, labelling from top left to bottom right. The lines represent the front, middle 
back edges of the WKB wave solution, where the pulse enters from the top of the box. 



(c) 
and 



dvy 
dx 



= 0, 



x=0.5 



dVy 



= 



-2.5 



dVy 



= 0, 



z=2.5 



Tests show that the central behaviour is unaffected by these 
choices. The other boundary conditions follow from the re- 
maining equations and the solenodial condition, V • Bi =0. 
We found that the linear Alfven wave travels in from the side 
boundary and begins to spread out, following the field lines. As 
the wave approaches the (left) separatix (that passes through 
X = -1), the wave thins (but keeps its original amplitude). The 
wave eventually accumulates very near this separatrix. This can 
be seen in Figure [H 

4.1.4. analytical results 

By solving the same characteristic equations, equations [141 but 
with diff'erent choices of initial condition (i.e. modelling the 
wave coming in from the side), we can obtain an analytical re- 
sult to match our numerical simulation. Again, we use a fourth- 
order Runge-Kutta method. The comparison of the numerical 
simulation and our WKB solution can be seen in Figure [H 

4.2. Two null points not connected by a separator 

We now consider the two null point magnetic configuration that 
does not contain a separator, as shown in Figure [T] (right). 



4.2.1. top boundary 

We consider a box (-2 < x < 2, -2 < z < 2) with a single wave 
pulse coming in across part of the top boundary (-1 < x < 1, 



z = 2), i.e. crossing the separatrices. The final, high resolution 
run was performed in a box (-2.5 < x < 2.5, -3 < z < 2). The 
full boundary conditions were; 



Vy(X, 2) ■■ 
dVy 

dx 



2(sina;0[sin^^] for 




0, 

OX 



-2.5 



-1 < JC< 1 

Otherwise 

dVy 

dz 



0, ^ 



. 



The other boundary conditions are chosen in the same way as 
before. 

We found that the linear Alfven wave travels down from the top 
boundary and begins to spread out, following the field lines. 
The wave thins as it descends, but keeps its original amplitude. 
The wave takes the shape of the separatrices that pass through 
X = ±2, z = + V2/3. The wave accumulates along these sep- 
aratrices for z > 0. The rest of the wave (the part between 
-I < X < continues descending and eventually accumu- 
lates along the separatrices in the lower half of the box. This 
can be seen in Figure [9l The part of wave in the small region 
about X = descends indefinitely, since the separatrices con- 
verge along X = at infinity. This can be shown: 



d 

'Jz 



Let ^ = {Bz§i) and comparing the original equation with = 



+ pleads to: 

OS ox OS Oz 



z = tan (A -01 



x=0 
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Fig. 8. Comparison of numerical simulation (shaded) and analytical solution at times {a) ^=0.14, {b) r=0.25, (c) t =0.5, {d) t=0J5, 
(e) t=l.25 and (/) ^=1.75, labelling from top left to bottom right. The lines represent the front, middle and back edges of the 
wave, where the pulse enters from the side of the box 



where A = arctan zo and zo is the starting postion of our char- 
acteristic. Thus z starts at x = 0, z = 2 and descends, following 
the line x = but becoming closer to zero. The particle passes 
through z = Sitt = arctanl and becomes more and more neg- 
ative (i.e. continues to descend) thereafter. 

4.2.2. analytical results 

Again, we have to solve equations ([T3l) to gain an analytical so- 
lution for our Alfven wave. Here, - 2xz and 8^ = x^-z^-l. 
We substitute V3; = ^^V'U,^) into ([13]) and assume that a; » 1 dvy 
(WKB approximation). This leads to a first order PDE of the dx 



form ( 



dz 



^ = 0. Applying the method of character- 



istics to this magnetic setup, we generate the equations: 



J0 

ds 
dp 

ds 
dq 

ds 
dx 

ds 
dz 
ds 



- (2zp + 2xq) f 

- (2xp - 2zq) f 
Oxz)^ 



(15) 



where f = ^xzp + (x^ - - l) = ^ = ^, 61; is the 
frequency of our wave and s is some parameter along the char- 
acteristic. 

These are solved using a fourth-order Runge-Kutta method and 
the results can be seen in Figure |9j 



4.2.3. side boundary 

We now consider a box (-2 < x < 2, -2 < z < 2) with a 
single wave pulse coming in across the side boundary (x = -2, 
< z < 1). The final, high resolution run was performed in 
a box (-2 < jc < 0.5, -2.5 < z < 2.5). The full boundary 
conditions were; 



Vy(-2, Z) 




-1 <Z< 1 

0<t<^ 

CO 

otherwise 



dVy 



0. 



z^2.5 



z^-2.5 

We found that the linear Alfven wave travels in from the side 
boundary and begins to spread out, following the field lines. As 
the wave approaches the separatix (that passes through x = - 1 , 
z = V2/3), the wave thins (but keeps its original amplitude). 
The wave eventually accumulates very near this separatrix. 
This can be seen in Figure [TOl 

4.2.4. analytical results 

We use the same characteristic equations, (equations [15]), with 
diff'erent choices of initial condition (i.e. modelling the wave 
from the side boundary this time) to obtain an analytical result 
to compare our numerical simulation. The results can be seen 
in Figure [To! 

5. Conclusions 

This paper describes an investigation into the nature of MHD 
waves in the neighbourhood of two null points. From the work 
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explained above, it has been seen that when a fast magne- 
toacoustic wave propagates near the two null point configu- 
ration, the wave bends due to refraction and begins to wrap 
itself around the null points. The wave 'breaks' and part of it 
travels (and wraps) around one of the null points and the rest 
wraps around the other null point. In the case of the fast wave 
approaching the two null points from above, the wave travels 
down towards the null points and begins to wrap around them. 
In addition, since the Alfven speed is non-zero between the null 
points, the wave can travel through this area (into the lower half 
plane). This part of the wave is also aff'ected by the refraction 
eff'ect and the wave continues to wrap around its closest null 
point. The wave breaks into two along the line x = (due to 
symmetry), with each half of the wave going to its closest null 
point. Each part of the wave then continues to wrap around 
its respective null point repeatedly, eventually accumulating at 
that point. In the case of the fast wave pulse travelling in from 
the side boundary, we see a similar eff'ect (i.e. a refraction ef- 
fect, wave breakage and accumulation at the nulls), but in this 
case the wave is not equally shared between the null points. 
For a fast wave travelling in from the left boundary, initially 
the pulse thins and begins to feel the eff'ect of the left hand- 
side null point. The wave begins to wrap around this null point 
(due to refraction). As the ends of the wave wrap around behind 
the left null point, they then become influenced by the second, 
right hand side null point. These arms of the wave then pro- 
ceed to wrap around the right null point, flattening the wave. 
Furthermore, the two parts of the wave now travelling through 
the area between the null points have non-zero Alfven speed, 
and so pop through. These parts of the wave break along x = 
and then proceed to wrap around the null point closest to them. 



So it is clear the refraction efl'ect focusses all the energy of 
the incident wave towards the null points. The physical signif- 
icance of this is that any fast magnetoacoustic disturbance in 
the neighbourhood of a null point pair will be drawn towards 
the regions of zero magnetic field strength and focus all of its 
energy at that point. 

From this work, it has been seen that when a linear, fast mag- 
netoacoustic wave propagates near a two null point configu- 
ration, the wave bends due to refraction and begins to wrap 
itself around the null points (at least in two dimensions). The 
angle that the fast wave approaches the null points from will 
determine what proportion of wave ends up at each null point 
(when the wave 'breaks'). Furthermore, since we have a chang- 
ing perturbed magnetic field with increasing gradients, we will 
have a build up of current density. In the case of the fast 
wave, all the wave is accumulating at the null points. This 
means that the perturbed magnetic field (Bi) will have large 
gradients at those points and that is where current will ac- 
cumulate. This is in good agreeement with similar phenom- 



ena noted in [McLaughlin & Hood (2004). With a large current 
accumulation at the nulls, this is where energy will be dissi- 
pated. Therefore, wave heating will naturally occur at these null 
points. Further experiments are being carried out to confirm 
this. 

In the cases of the Alfven wave, the results show that the wave 
propagates along the field lines, thins but keeps its original 
amplitude, and eventually accumulates very near to the sepa- 
ratrices. This is again in agreement with work carried out for 
a single null by McLaughlin & Hood (2004 ), As hypothesised 
above, since we have wave accumulation along the separatri- 
ces, the gradients of the perturbed magnetic field (jx and 
in this case) will build up along these as well. Hence we will 
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Fig. 10. Comparison of numerical simulation (shaded area) and analytical solution at times {d) r=0.14, (h) ^=0.39, (c) r=0.75, (d) 
t-XS), {e) t=l.25 and (/) ^=1.75, labelling from top left to bottom right. The lines represent the front, middle and back edges of 
the wave, where the pulse enters from the side of the box. 



have current build up along the separatrices. Therefore, if all 
the current is accumulating along the separatrices, then this is 
where dissipation will naturally occur. Our early experiments 
into current accumulation for two null points confirm this. In 
view of the fact that the waves eventually dissipate due to re- 
sistive damping (for example), it would be useful to have an 
estimate of the nature of the wave damping. Dissipation enters 
our model through the addition of a t/V^Bi term to equation d?]) 
and the WKB model can be extended to include the eff'ects of 
resistivity. For the model describing the Alfven wave pulse en- 
tering the top boundary (Section B.l.lb , it was found that along 
the separatrices the diff'usion term become important in a time 
that depends on log 7/ (see Appendix); as found by Craig & 
McClymont (1991, 1993) and |Hassam (T992), This means that 
the linear wave dissipation will be very efficient. This method 
can also be applied to the other cases investigated in this paper. 

Demoulin et al (1994b) investigated the magnetic field topol- 
ogy of a flare event. They found that the high intensity regions 
of Ha were located on or close to the separator lines. Their 
findings show a link to ours and support the possibility of trans- 
ferring our results to 3D. We can also easily extend the model 
to 2.5D with the addition of a third spatial coordinate; by tak- 
ing into account an extra Fourier component of the form e^^^, 
where m is the azimuthal mode number. This would lead to 
coupling of all the wave modes, and probably result in energy 
accumulating at both the separatrices and the null points. 

Note that all the experiments conducted here that send in an 
Alfven wave pulse from the boundary do so across a separa- 
trix. If the initial wave pulse does not cross a separatrix, it will 
follow the field lines and leave the box (there will be no wave 
accumulation at any separatrices). 



Appendix 




0.0 0,1 0.2 0.3 0.4 0.5 

s 

Fig. 11. Behaviour of log ( 00 against time elapsed. The slope 
of the line between ^ = 0.1 and s = 0.5 is 4oj. 



We can learn something about the nature of the wave damp- 
ing by considering an extension to our WKB expansion for 
the Alfven wave investigated in Section (14.1.1b . We assume 
solutions of the form e^^^'^~^^\ where » 1 and is (now) 
of order unity. We expand in terms of l/oj, such that cp = 
00 + ^01 + ^ 02 + ••• and consider r] = t]q/oj^, where tjq is 
also of order unity. Here we have assumed a specific scaling 
for the resistivity. This is for illustration only; alternative forms 
are possible and they modify the series expansion for 0. 
Using the same method with which we reached equation (fTTI) . 
the linearised equations for the Alfven wave including resistiv- 
ity can be combined to form a single wave equation for by\ 

d^by r. rjo d o 
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Substituting in by 



-o/by = -o/(Bo-V(f)fby-\-iojby(Bo'Vf(f) 
+ 7]o{ioj\V(pfby-\-byV^(p) . 

Now let = 00 + ^^01' where we assume 0i has a complex 
form (for consistency later) and we have taken out a factor of 
T]o (for simplification). 

-oj^ = -a? (Bo • V0o)^ - 2icom (Bo • V0o) (Bo • V0i) 
+ 7]l (Bo • V0i)2 + ioj (Bo • 00 - m (Bo • V)2 01 
+ /77o6i;|V0oP + [other terms] . 

Comparing terms of order of co^ and co: 

(Bo • V0o)' = 1 ^ (Bo • V0o) = 1 , 
-27/0 (Bo • V0o) (Bo • V0i) + (Bo • V)' 0o + m |V0ol' = 



•(Bo-V)0i 



2 |V0or 



Now: 



/^^o ^ ^00 



500 \^ /50O 



where P = ^ and ^ = ^ as in section (14.1.21) . Hence: 



dx 

J01 501 501 ^ f 2 2\ 

//v //v fir //v /97 9 ^ ^ ^ 



dz 
dz 501 



J^s" dx ds dz 



(16) 



where ^ = [(x^ - - l) - 2xz^], and we have used the 
characteristic equations from section (14.1.21) . Hence, we can 
solve for 0i by adding equation ([T6l) to ([T4l) and utilising the 
same Runge-Kutta method. The behaviour of 0i can be seen in 
Figure [TT] where log0i along the separatrix is plotted against 
s. The behaviour is described by 0i = ^e^^^^ where co - In. 
Hence: 

by = exp [/ (66)00 - ^0] exp [-(jj^ri^ij , 

= exp [/ (66)00 - ^0] exp |-^6l)^7/ [exp (46l)^)]| . (17) 

Therefore, the diff'usive term will become important when 
^(j?rie'^^^ becomes of order unity, i.e. on a timescale of ^ ^ 
-j^ log 7/, under the above assumptions. 
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